Phonon Dispersion Effects and the Thermal Conductivity Reduction in GaAs/AlAs 

Superlattices 



O 

o 



W. E. Bies, R. J. Radtke,* and H. Ehrenreich 

Physics Department and Division of Engineering and Applied Sciences, 
Harvard University, Cambridge, Massachusetts 02138 
(February 1, 2008) 

The experimentally observed order-of-magnitude reduction in the thermal conductivity along the 
growth axis of (GaAs)„/(AlAs) n (or nxn) superlattices is investigated theoretically for (2x2), (3x3) 
and (6x6) structures using an accurate model of the lattice dynamics. The modification of the phonon 
dispersion relation due to the superlattice geometry leads to flattening of the phonon branches and 
hence to lower phonon velocities. This effect is shown to account for a factor-of-three reduction 
in the thermal conductivity with respect to bulk GaAs along the growth direction; the remainder 
is attributable to a reduction in the phonon lifetime. The dispersion-related reduction is relatively 
insensitive to temperature (100 < T < 300K) and n. The phonon lifetime reduction is largest for the 
2x2 structures and consistent with greater interface scattering. The thermal conductivity reduction 
is shown to be appreciably more sensitive to GaAs/AlAs force constant differences than to those 
associated with molecular masses. 
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PACS numbers: 63.22.+m, 84.60. Rb, 66.70.+f, 81.05.Ea 



I. INTRODUCTION 



Superlattice. structures haverJpeen proposed to be materials with a high thermoelectric figure of merit ZT, for 
both in-planca and cross-planecra current flow. In the latter case, the improvement in thermoelectric performance 
is attributable to a reduced lattice thermal conductivitytTta rather than to a higher electronic conductivity. Exper- 
imentally, a factor-of-ten reduction in the component of lattice thermal conductivity along the growth axis, Ke.zz, 
is observed in GaAs/AlAsBB and Bi2Te3/Sb2Te3L]tl superlattices (iSLs). Theoretically, the thermal conductivity of 
Si/Ge SLs has been studied previously by Hyldgaard and MahanO and by ChenE2l. Within the context of a very 
simplified model of the phonon dynamics, the calculations of Ref. 9 were able-4o reproduce the factor-of-ten reduction 
in the thermal conductivity along the growth direction. By contrast, Chen'sEd extensive work focussed on the role of 
thermal boundary resistance. The present paper investigates the reduction associated with SL induced changes of the 
I/"") . phonon dispersion based on a realistic, computationally intensive treatment of the phonon spectra and dynamics. 

The origin of the observed reduction in thermal conductivity may be explained qualitatively as follows. First, we 
note that, in the experimental work of Capinski and MarisO and Capinski et alu on (GaAs)„/(AlAs) n SLs for n up 
to 40, the phonon mean free path inferred from the thermal conductivity, heat capacity and Debye velocity is greater 
than 370 A at all termperatures for which measurements exist, always large compared to the size of the SL unit 
cell, 5.66 A. Thus, the phonon transport lies in a regime where the SL phonon dispersion relation and lifetime, and 
not those of the bulk constituents, determine the thermal conductivity. According to the expression for the thermal 
conductivity derived from the phonon Boltzmann equation in the relaxation-time approximation (see Eq. (6)), the 
thermal conductivity depends on (1) a quantity representing the contribution of the SL phonon dispersion relation, 
and (2) the lifetime r, which contains phonon-phonon, interface and defect scattering effects. We shall focus only on 
O ■ item (1), whose effect can be computed within our realistic, albeit complicated, lattice-dynamical model. The effect 
, of the SL geometry is to introduce anticrossings and new gaps in the phonon dispersion relation, when the magnitude 
of the phonon wavevector along the growth direction equals an integer multiple of n/d, where d is the period of the 
SL along the growth direction. The consequent flattening of the phonon branches near the Brillouin zone edge leads 
to a lowering of phonon velocities in the growth direction, and hence a reduction in thermal conductivity. 

We describe in Sect fop, IL the lattice-dynamical model for the SL, which is a generalization of the 11-parameter 
rigid-ion model of KuncPH 12 !. It incorporates short-range interactions to next nearest neighbors and the long-range 
Coulomb interaction. The construction of the SL dynamical matrix is outlined in Section II and in the Appendix. 

The formalism is applied to a (GaAs)3/(AlAs)3 SL in Section III. The phonon dispersion relation will be seen to 
display the flattening expected for a SL. Critical points, especially at the high-symmetry points T, X and Z, produce 
sharp peaks in the density of states in the SL. Miniband formation and anticrossings in the SL phonon dispersion 
relation lead to a three-fold reduction, relative to bulk, in the contribution of the phonon dispersion relation to t 
thermal conductivity along the growth direction. The present results contrast with those of Hyldgaard and Maha: 
who found that, in their simplified model, the order-of-magnitude reduction of k zz /t was attributable to effects 
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related to the SL dispersion relation alone. In our more realistic treatment of the SL dynamics, a significant three-fold 
reduction in the lifetime is also needed to explain the experimental reduction in n zz by a factor of ten. Finally, the 
sensitivity of the decrease in n zz to differences in masses or force constants between the GaAs and AlAs layers is 
investigated; differences in the force constants are found to play a markedly greater role in the reduction of k zz than 
differences in mass. 

Section IV is devoted to discussion and conclusions which are broadened by examining the dependence of our results 
on the period n of the (GaAs) n /(AlAs)„ SL for n = 2, 3 and 6. These results permit some discussion of interface 
effects on the phonon lifetime and a more detailed comparison of the present work with that of Refs. 9 and 10. The 
reduction in the contribution of the phonon dispersion relation to the thermal conductivity of the SL relative to bulk 
is computed, and found to be approximately independent of n for the small values considered here. The lifetime for 
both the SL aad bulk was determined using the experimental thermal conductivities of the SL and of bulk GaAs of 
Capinski et al H. The lifetimes are found to be smaller for the 2x2 SL and larger and roughly equal for the 3x3 and 
6x6 SLs, consistent with the presence of greater interface scattering in the 2x2 SLs. 



II. FORMALISM 

The lattice dynamics can be treated realistically via the 11-parameter rigid- ion model of KuncSIll, which has been 
successfully applied to zinc-blende-structure compounds. Consider (GaAs)3/(AlAs)3. In the SL unit cell, consisting 
of a layer of GaAs above a layer of AlAs in the growth direction, there will be 3 pairs of GaAs followed by 3 pairs 
of AlAs, or 3 Ga, 3 Al and 6 As atoms in all, which may be indexed by K = 1, . . . , 12. Letting denote the 

displacement in the direction a — x, y, z of the K-th atom in the ^-th unit cell, plane-wave solutions of the form 

Ua (£ K ) = M- 1 / 2 e i(k ' x( ^ ) -^ (k)t) u; a (K|ki) (1) 

are assumed, where x(£k) is the equilibrium position of the K-th atom in the ^-th unit cell, and w a (K\k.j) satisfies the 
secular equation 

^■(k) 2 WQ ( K |kj) = ^2c a p(KK'\k)w p {^\kj). (2) 

/3k' 

The dynamical matrix C_ reflects the interatomic force constants of the crystal. In the present rigid-ion model, the 
interatomic forces are divided into (1) short-range forces extending to second nearest neighbors, and (2) the long-range 
Coulomb interaction. Accordingly, the dynamical matrix may be written 

Q. = Q sr + £Zcoul- (3) 

As a result of symmetry of the zinc-blende structure, the short-range forces to second nearest neighbors may be 
described by ten parameters for each materialE3. For nearest and next nearest neighbor interactions in the SL 
unit cell, we employ the force constants determined for the constituent bulk materials separately,, rot ated by the 
appropriate point-group operation. Bulk GaAs and AlAs parameters are taken from the literaturetSlla. In the SL, 
bulk parameters are used within each layer. For the interface atoms and Ga-Al bonds crossing the interface, we 
employ the average of the bulk parameters following Ren et at-3. For the Coulomb interaction, the atoms are treated 
as point charges. The Madelung sum, and its derivatives, are computed using the usual Ewald transformation!!^, 
which has been generalized here for SLs. This is accomplished by separating the sum over the spatial index £ into a 
sum over layers normal to the growth direction, and a sum along the growth axis, £ z . A two-dimensional Ewald 
transformation is performed in each layer; these results are then summed over £ z . The resulting expressions for the 
function 0(k, r) and its derivatives (see the Appendix) are similar to those that arise in the usual three-dimensional 
Ewald procedure; however, the definite integrals differ and must be performed numerically. Detailed expressions for 
the Coulomb term are given in the Appendix. 

The phonon Boltzmann equation in the relaxation-time approximation leads to the following expression for the 
lattice thermal conductivity: 

d3q ^^^^^T ph (<>,T) (4) 



(2tt) 3 ^ q % % dT 

where n(u^) is the distribution function of the phonons, the sum is over branches a, and r p h(wq a \ T) is the lifetime. 
Eq. (4) can be written in terms of 
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(5) 



dqi dqj 



as 
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dw(dn(o;)/dT)Sy(a;)rph(w, T). 



(6) 



This paper will focus on the SL effects on the dispersion relation contained in E^ (u>) . Relaxation-time effects associated 
for example with scattering from interfaces, defects, umklapp processes, etc. are not considered explicitly. However, 
the results will be used to infer some of their properties. 



We focus first on the (GaAs)3/(AlAs)3 SL studied experimentally by Capinski and MarisQ. Using a picosecond 
pump-and-probe technique, they observed an order-of-magnitude reduction in the thermal conductivity along the 
growth direction, k zz , relative to bulk GaAs. The dispersion relation along the TX and TZ directions, which was 
generated numerically according to the method described in Section II, is shown in Fig. 1. The significance of the 
labelled features will be explained in the discussion of Fig. 2. Because the SL unit cell contains three unit cells each 
of bulk GaAs and bulk AlAs, respectively, arranged along the growth axis, the edge of the SL Brillouin zone in the 
growth direction, Z, is one-sixth as far from the center as it is in the in-plane directions, X and Y. As a result each of 
the six branches in the bulk material is folded back six times along TZ, which is most easily seen for the longitudinal 
acoustic mode in Fig. 1. Bulk GaAs and bulk AlAs optical modes have no frequencies in common, and are therefore 
localized. This leads to (1) flat SL dispersion in optical modes, and (2) localization of AlAs optical modes to the AlAs 
layer. The GaAs optical modes are not localized, since they overlap with the acoustic modes as shown in Fig. 1 along 
TX. Due to the non-analyticity of C Coul as q — > in the SL, uj(\q\ — > 0) differ along TX and TZ (See Ref. 14). 

The density of states p{oj) versus frequency, computed using the tetrahedral integration method as presented in 
MacDonald et al liZI, is given in Fig. 2. Note the transverse-acoustic features around ui = 80 to 100 cm -1 , the 
longitudinal-acoustic features around 150 to 200 cm -1 , the GaAs optical feature at 220 to 260 cm' 1 , and, separated 
in frequency at higher frequencies, the AlAs optical feature at 330 to 400 cm -1 . Band gaps at the zone edge and 
anticrossings in Fig. 1 yield critical points which, depending on the amount of q-space at those frequencies, produce 
sharp structure in the density of states. As shown in Fig. 2, this structure in the density of states can be correlated 
with features in the dispersion relation, usually at the Y (for folded-back bands), X and Z points. The structures at 
r, X and Z labelled in Fig. 2 are identified with the corresponding features in Fig. 1. The strength of each feature 
depends on the integral of 5(u> — Wq^) (cf. Fig. 2) at that frequency, which will be large if |v| is small. Surprisingly, 
most of the peaks in the density of states appear to be associated with the critical points at Y, X and Z. No such-fine 
structure in the density of states exists in bulk GaAs or bulk AlAs, as may be seen for instance in Patel et al EJ for 
GaAs. (The density of states for AlAs is not available in the literature, but our calculations confirmed the absence of 
fine structure for AlAs as well.) 

Fig. 3 shows the results for E(uj), as defined in Eq. (5), for bulk GaAs and the SL. Note that optical modes do 
not contribute appreciably to T, zz (ui) in the SL, an effect of localization in the AlAs layers: flat dispersion leads to a 
vanishing of dto^ /dq z . The fine structure in the density of states is also correlated with that in E(w): peaks in the 
density of states DOS oc J dS/\v\ become dips in E cx J v 2 dS/\v\. (Here, S denotes the surface of constant frequency 
lu in the Brillouin zone; it consists not only of the surface containing the critical point, but also possibly of other 
surfaces at the same u> elsewhere in the zone.) Acoustic modes in the SL contribute less than they would in bulk 
because of the band-gap and anticrossing-induced reduction in v 2 . This leads to a three-fold reduction in ki/t at 
300 K, determined here either by integrating Eq. (4) directly on a 60 x 60 x 20 grid oavering an irreducible wedge 
of the Brillouin zone, or by integrating Eq. (6). Experimentally, Capinski and Marisu found a ten-fold reduction 
factor for (GaAs)3/(AlAs)3. The full reduction factor is a product of the reduction due to E and that due to r. 
Assuming r to be constant at any given temperature, it can be found by requiring equality in Eq. (4) or (6) to the 
experimental value for the thermal conductivity in bulk or SL, respectively, if the appropriate dispersion relation is 
used in computing k zz /t. The lifetime in bulk versus lifetime in SL, as well as the layer-width dependence of the 
results, will be discussed below. Finally, we note that T, zz < Yi xx in the SL for the simple physical reason that band 
flattening along the ^-direction affects du/dq z more than du/dq x . 

The reduction in thermal conductivity due to E can be understood by means of an easily visualized picture in 
q-space. In bulk, the quantity q x ^ a {dn{io^) / dT)io^\v^],) 2 represents the contribution of the phonon dispersion 
relation in Eq. (4). In a range Aq x of the integrand it is weighted by q x because, the dispersion relation being 
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rotationally symmetric in the (q x ,q y ) plane, we may integrate around the circle of radius q x to yield a properly 
weighted function of q x and q z alone. The quantity is plotted in Fig. 4(a),(b),(c) (bold line) together with the 
corresponding values for GaAs (light solid line) as a function of q z for three values of q. x . The SL contribution is 
reasonably localized in q z . To gain physical insight, we replace it by the dashed rectangles (of equal area). This 
approximate localization is related to the reduction at q z rs and ir/d due to miniband formation (that is, flattening 
of lu vs. q) and has dips from anticrossings, as in Fig. 4(a), for instance. The dependence on q x can then be 
summarized by the density plot in Fig. 4(d), comparing the SL on the left to bulk GaAs on the right. For GaAs 
the equivalent rectangles extend over the entire range of q z because there is no localization due to band flattening at 
the zone edge as in the case of the SL. The shading indicates the weight of each increment Aq x . We find that points 
around q x = w/d and at the zone edge (q x = dir/d) contribute most to heat transport. The latter is due to the effect 
of the weighting by q x in the annular integration. In addition, this weighting causes the contribution from points near 
the origin to be very small. A simple numerical estimate leads to a reduction in E to 34%, which is to be compared 
to 36% in the exact calculation. 

Finally, we discuss the sensitivity of variation of the mass differences, force constants and effective ionic charge 
e* between layers on the thermal conductivity. This manifests itself through variations of the zone edge gap and 
hence the group velocity. This effect has been studied by interpolating e* , the force constants K, and cation mass M 
between their AlAs and their GaAs values in the GaAs layer. Here K refers collectively to the 10 Kunc parameters 
describing interatomic forces. Explicitly, starting from (AlAs)6 we (i) vary e* 2 and M linearly in three neighboring 
AlAs positions to make a (AlAs) 3 /(GaAs)3 SL. In the case of the masses we put M = (1 — a)Mx\ + cxMa a so that 
at a — and 1 correspond respectively to Al and Ga; (ii) vary e* 2 and force constants K linearly in a in the same 
way; and (iii) vary e* 2 , and both M and K linearly in a. The results are given in Fig. 5. The thermal conductivity 
is seen to be far more sensitive to the variation of force constants than the variation in mass. The variation of the 
force constants alone produces band flattening which reduces the thermal conductivity by about 60%. The additional 
flattening due to changing masses leads to only a few percent additional reduction. 



IV. DISCUSSION AND CONCLUSIONS 



Before turning to the final results, we emphasize again that this paper is primarily concerned with the effects of 
superlattice induced changes of the phonon dispersion on the lattice thermal conductivity. Since a SL is a perfect 
crystal, the ideal structures under discussion here automatically include the coherent effects associated with perfect 
interfaces. The importance of including such effects was first pointed out by,Hyldgaard and Mahand in connection 
with a simple, highly idealized model for Si/Ge SLs. Extensive work by ChenEj focussed on diffuse interface scattering 
of phonons. His model assumes SL layers sufficiently thick that the phonon spectrum in each layer corresponds to 
that of the bulk. A mixture of spectral and diffuse interface processes is found to be sufficient to explain the observed 
experimental reduction of the SL thermal conductivity relative to bulk. As already pointed out, the more realistically 
modeled results of the present paper for Kg/r lead to a three fold reduction without any assumptions abcoit the SL 
phonon scattering mechanisms. Our results must therefore be viewed as complementary to those of ChenE3. Taken 
together, they suggest that a combination of phonon spectral changes and imperfect interfaces can account adequately 
for the observed reduction. 

The present calculations permit some statements concerning lifetime effects from the dependence of the SL results 
on layer width. This dependence was studied numerically for 2x2, 3x3 and 6x6 GaAs/ AlAs SLs. The results for 

A y- . ( Q ) du4 a) &4 a) dn(J™ ] ) 



(2tt)3^ q dq t dqj dT 

are given in Table I. r is assumed constant. Given an experimental value for n zz and the calculated value of E Z2 , the 
lifetimes listed in Table I are found as r = k zz /T, zz ; only the lifetime itself is given in Table I and not E 22 , but the 
ratios E 22 (SL)/E Z2 (bulk) are listed because we are interested in 

r SL _ E 22 (bulk) >4 L (expt) 



Tbuik E 22 (SL) K^ lk (expt)' 

The SL E 22 is found to have a value about 40% of the bulk, and to be relatively insensitive to the SL period and 
temperature. For larger SL periods, there is more folding back, but this is compensated by the smaller size of the 
gaps at the zone edge, which is found in the computed TZ dispersion relations to scale inversely with the SL period. 
The two effects balance, resulting in an approximately constant reduction factor. The phonon lifetimes for the 2x2, 
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3x3 and 6x6 SLs are also given in Table I. The bulk lifetime is the same to within 2% for the different SL periods, as 
it must be. The ratio of SL to bulk lifetimes is significantly smaller for the 2x2 SLs at each temperature while it is 
larger and roughly the same for the 3x3 and 6x6 SLs. Thus, the reduction in thermal conductivity may be divided 
into a dispersive part which is insensitive tonxn and a scattering part which is sensitive to interface scattering. The 
results for the 2x2 SL are possibly associated with the experimental difficulties associated with achieving sufficiently 
perfect interfaces for small n. 

The present calculations imply: (1) that a similar reduction in the contribution of the SL phonon disperion relation 
to transport in the growth direction, and perhaps in k zz itself, may be expected in any SL with similar mass or 
force-constant differences between layers. (2) If the lifetime is reduced in SLs by increased umklapp scattering, as 
suggested by Ren and Dow 13 , then the present calculations give an upper bound on the SL n zz and a lower bound 
on the increase in ZT in a SL relative to bulk. 
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APPENDIX A: 



This Appendix presents the detailed formulae for the Coulomb part of the dynamical matrix, derived using the 
Ewald procedure as described in the text of the paper. Letting k label the atoms of the SL unit cell, M K be the mass 
of the K-th atom and x(kk') be the separation vector from the K-th. atom to the n'-th atom, we have, for k ^ n' , 



aP V 1 ' dr a dr 



<.(k'k) 



where 0(k, r) is by definition 



0(k,r)=£ 



|x(*) + r| 



(Al) 



(A2) 



x(£) being the position of the ^-th unit cell. The result of the Ewald procedure adapted to the superlattice is that 
4>(k, r) can be written in the form 

0(k, r) = R^2H(\x(£) + r|i?)e* k - x W 



2^ 



/(0O) W+M Kftn) + k nii x (^) + A )c -i(T(h„Hk,|Vr„ + i fc ,s.»(^ (A3) 



+ ?-f |r(/i||)+k||r"' 2R 2 



where R is an arbitrary cutoff (we always take R = 3/oo), 



H(x) 



I 

J X 



(A4) 



uii is the area of the unit cell in the superlattice plane, £ z labels the layers perpendicular to the growth (z) axis, 
= (h x ,h y ) labels the cells located at reciprocal lattice vectors r(h\\) in each plane, k = (ky, k z ), and 



7(a,/3,7) = / e- v2 -* 2 / v2 dv. 
J0 



For k — k' we have 



ik-x(£) 



dr a drp 



dr a drp 



t=x(£'k',0k) 



(A5) 



(A6) 
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The first term on the right is similar to the above expressions in Eqs. (A1)-(A3) but for the I = in Eq. (A3) term 
one substitutes (4/3y/n)S a fj. The second term is given by 

<h [/0+J1+/2U0 (A7) 



dr a drp 
with 

/„ = Re K H°(\r\R), (A8) 



h=R e K ,H(\x(e')+x( K 'K) + r\R) (A9) 



and 



where 



A similar, but not identical, approach was used in Ref. 15. 



H o {x) = _vvi 1 e - x - dx , (A11) 



* Present address: Schlumberger, Sugar Land Product Center, 110 Schlumberger Dr., MD110-4, Sugar Land, TX 77478. 
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TABLE I. Reduction in E zz = k zz /t relative to bulk GaAs for 2x2, 3x3 and 6x6 SLs at T = 100, 200 and 300 K, and SL 
phonon lifetimes deduced from Eq. (8) for the 2x2, 3x3 and 6x6 SLs at each temperature. 





experimental 


theoretical 








G T 


Kzz (W/cm KJ 


S«(SL) 
£z Z (bulk) 


Tbulk (PSJ 


TSL (PS) 


TSL/Tbulk 


T=300 K, K^" lk =0.45 W/cm K " 












2x2 


0.040 


0.38 


36.8 


8.7 


0.24 


3x3 


0.068 


0.36 


37.2 


15 


0.42 


6x6 


0.053 


0.34 


37.2 


13 


0.35 


T=200 K, K^ lk =0.64 W/cm K b 












2x2 


0.055 


0.38 


55.2 


12 


0.22 


3x3 


0.090 


0.39 


56.1 


20 


0.36 


6x6 


0.072 


0.35 


55.7 


18 


0.32 


T=100 K, K^" lk =2.0 W/cm K " 












2x2 


0.065 


0.40 


222 


18 


0.08 


3x3 


0.110 


0.42 


227 


30 


0.14 


6x6 


0.096 


0.37 


222 


29 


0.13 



a Experimental values for GaAs/AlAs reported by Ref. 
b Experimental value for GaAs listed in Ref. H. 
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FIGURES 



FIG. 1. (GaAs)3/(AlAs)3 SL dispersion relation along the rX=(27r/ao, 0, 0) and rZ=(0, 0, 7r/3ao) directions; a is 
the conventional GaAs unit cell size. The labels a — w are defined in Fig. 2. 

FIG. 2. Phonon density of states for the (GaAs)3/(AlAs)3 SL, with labelled critical points identified with features 
in the dispersion relation in Fig. 1. 

FIG. 3. The transport quantities Y* xx (oj) and T, zz (ui), defined by Eq. (5) in the text, for bulk GaAs (solid line) and 
the (GaAs) 3 /(AlAs) 3 SL (dashed and bold lines). 

FIG. 4. The transport quantity q x {dn/dT) J2 a i v q^z) 2 for < q z < ir/d at fixed q x , for (a) q x = 2n/d, (b) 
q x = 4w/d and (c) q x = 6w/d; solid line: bulk GaAs; bold line: (GaAs)3/(AlAs)3 SL. (d) Density plot in the (q x ,q z ) 
plane whose shading indicates the weight of each increment Aq x along q x to the value of this transport quantity for 
the (GaAs) 3 /(AlAs) 3 SL and bulk GaAs. 

FIG. 5. Dependence of n zz as e* 2 and the mass M (solid line), the spring constants K (long-dashed line) and both 
M and K (short-dashed line) are interpolated between their AlAs (a = 0) and GaAs values (a = 1) for the atoms in 
three contiguous layers in what is initially (AlAs) 6 at a = 0. 
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(GaAs) 3 /(AlAs) 3 SL 
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Fig. 1, Bies et ai, J. Appl. Phys. 



(GaAs) 3 /(AlAs) 3 SL 
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Fig. 2, Bies et ai, J. Appl. Phys. 
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Fig. 3, Bies et ai, J. Appl. Phys. 
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Fig. 4, Bies et ai, J. Appl. Phys. 
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Fig. 5, Bies et ai, J. Appl. Phys. 



